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Abstract. The notion of Bezout matrix is an essential tool in studying broad 
variety of subjects: zeroes of polynomials, stability of differential equations, 
rational transformations of algebraic curves, systems of commuting nonselfad- 
joint operators, boundaries of quadrature domains etc. We present a survey of 
several properties of Bezout matrices and their applications in all mentioned 
topics. We use the framework of Vandermonde vectors because such approach 
allows us to give new proofs of both classical and modern results and in many 
cases to obtain new explicit formulas. These explicit formulas can significantly 
simplify various computational problems and, in particular, make the research 
of algebraic curves and their applications easier. In addition we wrote a Maple 
software package, which computes all the formulas. For instance, as Bezout 
matrices are used in order to compute the image of a rational transformation 
of an algebraic curve, we used these results to study some connections be- 
tween small degree rational transformation of an algebraic curve and the braid 
monodromy of its image. 



Introduction 

Numerous works in operator theory shows that an algebraic curve given by a 
dcterminantal representation can be associated to a system of commuting non- 
sclfadjoint operators. In the simplest case there are two commuting nonsclfadjoint 
operators in a system and they are rational images of the same operator. This 
particular case gives a motivation and a framework for the studying the image of a 
complex line under a rational transformation. Using the notion of a determinantal 
representation of an algebraic curve there was found the explicit formula which 
describes the image of a complex line under a rational transformation. The image 
is given by a determinantal representation which uses the Bezout matrices of pairs 
of polynomials that define the rational transformation. 

This property of Bezout matrices can be used for many purposes, i.e. to study 
the braid monodromy of the image of two intersecting lines under rational transfor- 
mations of small degrees. Also using this property one can describe explicitly the 
boundary of a quadrature domain. This property and many others can be proved 
using the Vandermonde vectors, which are a natural framework to study Bezout 
matrices. 

These properties were used by the authors of this paper to create a package of 
procedures for Maple software This package allows the computation of all 
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classical and modern formulas that are mentioned in the paper, and to perform 
many different tasks. 

1. BEZOUT MATRICES AND VANDERMONDE VECTORS 

First example of Bezout matrices for polynomials of small degree appeared in 
Euler's work in 1748, jZj. Using this example Bezout gave a general definition 
of Bezout matrices for polynomials of any degree in 1764, j5|. The notation of 
Bezoutian matrix was introduced by Sylvester in 1853, |23j . The most common 
definition was given by Cayley in 1857, 0]. 

Definition 1.1. For two polynomials in one variable p(x) and q(x) of degree n 
there exists uniquely determined nXn symmetric matrix B(p,q) = (£>y)j -=i such 
that 

p(x)g(y) - q(x)p(y) 



x-y 



The matrix B(p, q) is called the Bezout matrix of polynomials p and q. 

We will use a framework of Vandermonde vectors. The matrices composed of 
Vandermonde vectors appeared inexplicitly in Vandermonde's work in 1772, |25j . 
The notation was attributed to Vandermonde by Weill in 1888, \27\ . 



Definition 1.2. 



V n {x) 



( 1 \ 

X 



V n (x) is called the Vandermonde vector of the length n. 

Theorem 1.3. If xi,X2, ■ ■ ■ , x n are pairwise distinct then vectors 
V n {xi), V n (x2), ■ ■ ■ , V n (x n ) are linearly independent. 

Proof. Let us consider n x n matrix 





I 1 1 . 


■ 1 \ 


(V n (xi), V n (x 2 ), ■ • ■ , V n (x n )) = 


X\ X2 


x n 




\ X-y X 2 


■ J 



Let us suppose that p = (po,pi, . . . ,p n ^i) is a row-vector from the left kernel of 
this matrix. Then the polynomial p(x) — po + pix + ■ ■ ■ +p„_ia; n_1 has n zeroes. 
The degree of this polynomial is less than n. Hence, the kernel of this matrix is 
trivial and vectors V n (xi), V n {x2), • ■ ■ , V n {x n ) are linearly independent. □ 

We will call V n {x) the Vandermonde vector of order zero. It is natural to de- 
fine the Vandermonde vector of higher orders. For higher orders the definition is: 

( x ) — yTEptZ 1 ) We will call V% (x) the Vandermonde vector of order k. 

Theorem 1.4. If ii + £2 + • • • +2 m = n — 1 and x\, X2, ■ ■ ■ , x m are pairwise distinct 
then vectors V n ( Xl ), Vfai), V£ (xi), V n (x 2 ), V^(x 2 ),..., V^(x 2 ), 
. . . , V n (x m ), V^(x m ), . . . , V^ m {x m ) are linearly independent. 



The prove is the same as of the above theorem. 
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Lemma 1.5. For every two polynomials p{x) and q{x) 
1. 



2. 
3. 



V 1 (x)B(p,q)V(y) 

x-y 

V T (x)B(p, q)V{y) = V T (x)(B(p, l)q(y) - B(q, l)p(y))V{y) 



V T (x)B(p, q)V(x) = q(x)p'(x) ~ p(x)q'(x) 

Proof. The first statement follows immediately from the definition of the Bezout 
matrix. The second statement follows from the next decomposition: 

p(x)g(y) - q(x)p(y) = p{x)g(y) - p{y)g(y) +p(y)g(y) - g(x)p(y) = 
x — y x — y 

t , p(x) -p(y) , N q(x) - q(y) 

q{y) p(y) = 

x — y x — y 

q(y)V T (x)B(p, l)V(y) - p(y)V T (x)B(q, l)V(y) 
To prove the third statement let us notice that 

V T {x)B(j>, l)V(x) = lim V T (x)B{p, l)V{x + e) = 

e— >0 

p(x)~p(x + e) p(x + e)-p(x) , 
lim ; : — = hm = p (x) 

e^O x - (x + e) e^O e 

This and the first statement of the lemma imply the third statement. □ 

Corollary 1.6. For every polynomials p(x) and q{x), vector w and point y 

w T B(p, q)V(y) = w T (B(p, l)g(y) - B(q, l)p{y))V{y) 

Proof. Vector w can be represented as a linear combination of Vandermonde vec- 
tors, w — ^2aiV(zi), where all points Zi are different from y. If follows from the 
second statement of the previous lemma that 

w Tnt„ - (sr. 



^B(p,q)V(y) = (J2^V T (z t )j B(p,q)V(y) = (V T (z t )B(p, q)V(y)) = 

2 at (V T ( Zl )(B(p, l)q(y) - B(q, l)p(y))V(y)) = 
ai V T ( Zi )) (B(p, l)q(y) - B(q, l)p(y))V(y) = 
w T (B(p,l)q(y) - B(q,l)p(y))V(y) 

□ 

One of the very well known classical applications of Bezout matrices is Jacobi- 
Darboux theorem. This theorem was proved independently by Jacobi in 1836, |11| 
and by Darboux in 1876, UJ. 

Theorem 1.7 (Jacobi-Darboux theorem). The number of common zeroes of two 
polynomials equals to the dimension of the kernel of the Bezout matrix of these 
polynomials. 

There are different proofs of this theorem. For the proof in terms of Sylvester 
matrices and Vandermonde vectors see [T%] . 
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2. Inverse of Bezout matrix and Hermite theorem 

In 1974 Lander proved that the inverse of Bezout matrix is a matrix of Hankel 
type, |14[ . This property, its applications and similar properties of structured 
matrices were studied by Gohberg and Olshevsky jS], Heinig and Hellinger [Sj, 
Tyrtyshnikov |24| . Chen and Yang 0. The framework of Vandcrmonde vectors 
allows not only to reprove the result of Landau but also to obtain explicit formulas 
for the coefficient of the inverse of Bezout matrix, • 

Theorem 2.1. If two polynomials p(x) and q(x) of degree n have no common zeroes 
and polynomial p(x) has no multiple zeroes then the inverse of Bezout matrix B(p, q) 
is a Hankel type matrix H — (/ly) and 

J+J-2 



k ^ q{xk)p' '(x k y 

where X\, x 2 ■ ■ ■ x n are zeroes of p(x). 

Proof. Let us construct a matrix V p from Vandermonde vectors in zeroes of the 
polynomial p(x): V p = (V(xi), V(x2), ■ ■ ■ , V(x n )). We will denote by D the ma- 
trix V p T B(p,q)V p . It is obvious that B(p,q) = (V p T y 1 D(V p )~ 1 and therefore 
B- 1 (p,q) = V p D- 1 V p T . 

Let us denote by dij coefficients of the matrix D. If i ^ j then 

d« = V T ( Xi )B(p,Q)V(xj) = P^i^-^M^) = 0> 

since the numerator is zero and the denominator is not because all zeroes of p{x) 
are different. If i — j then by the third statement of the Lemma H31 

<ki = q{xi)p'(xi) - p{x.j)q' \xi) = q(xi)p'(xi), 

therefore I? is a diagonal matrix: D = diag(q(xi)p'(xi), . . . , q(x n )p'(x n )) which 
means that 



D 1 = diag ( — — — 



\q(xi)p'(xi) q(x n )p'(x n ) 



Therefore 



V n- 1 = ( I 

[q((x k )p'(x k )). k 
Clearly, V p = {x]7 l )kj Hence, the inverse of Bezout matrix is 



/ n i+j-2 

B~ l {jp, q) = V p D- l V£ = I ' 



v , q{xk)p'{x k ) I 

which proves the theorem. □ 

Hermite theorem was proved by Hermite in 1856, jl()j . The theorem determines 
when all zeroes of a polynomial belong to the upper-half plane. This theorem can 
be proved in a framework of Vandermonde vectors, |19|. 

Theorem 2.2 (Hermite theorem). All zeroes of the polynomial p{x) belongs to the 
upper-half plane if and only if the matrix ^B(p,p) is positive definite. 
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3. Modern applications 



The simplest and very illustrative case of rational transformations of algebraic 
curves is a of a rational transformation of the complex line C into the complex 
plane C 2 . The image is a rational plane algebraic curve and the explicit formulas 
for this image in terms of the polynomials that define the rational transformation 
were obtained by Kravitsky in 1979, |13| . 

Theorem 3.1. Three polynomials in two variables po(x), Pi(x) and P2{x) map 
complex line C into complex plane C 2 : 



\Po(x)' p (x) J 

The image is a rational curve defined by a polynomial 

A(xi,X2) = det(B(pi,p 2 ) + xiB(p 2 ,po) + x 2 B(p ,pi)) 

Proof. Let x and y be two points on the curve, V{x) and V(y) be two Vandermonde 
vectors. V T {x)(p (x)B(p 1 ,p 2 ) + p 1 (x)B(p Q ,p 2 ) + p 2 {x)B(p ,pi))V(y) = 
p (x)V T (x)B(p 1 ,p 2 )V(y)+p 1 (x)V T (x)B(p 2 ,p Q )V(y)+p 2 (x)V T (x)B(p ()l p 1 )V(y) = 
(x - y)(po(x)(pi(x)p 2 (y) -P2(x)pi(y)) + (P2(x)p (y) - Po(x)p 2 (y)) + (po(x)p 1 (y) - 
Pi(x)po(y))) = This identity holds for for arbitrary y. Hence, (p (x)B(p 1 ,p 2 ) + 
pi(x)B(j)2,Po) + p 2 (x)B(p ,pi))V(x) = and therefore 



3.1. Nonselfadjoint operators. The work of M. S. Livsic and his collaborators 
in operator theory associates to a system of commuting nonselfadjoint operators an 
algebraic curve (called the discriminant curve) given by a detcrminantal represen- 
tation, see [TJI- This discovery leads to a very fruitful interplay between operator 
theory and algebraic geometry: problems of operator theory lead to problems of 
algebraic geometry and vice versa. 

A natural problem in operator theory is to define properly the notion of a ratio- 
nal transformation of a system of commuting nonselfadjoint operators. This arises 
whenever one wants to study the algebra generated by a given system of commut- 
ing nonselfadjoint operators. It may also allow representing the given system of 
commuting nonselfadjoint operators in terms of another system which is simpler 
in some sense (e.g., it contains fewer operators, or the operators have a smaller 
nonhcrmitian rank). A related problem in algebraic geometry is to find an im- 
age of an algebraic curve given by a determinantal representation under a rational 
transformation. 

To formulate the main problem more precisely, we have to introduce some nota- 
tion. We shall use the framework of commutative vessels ^3] which turns out to be 
very convenient in the study of commuting nonselfadjoint operators; it generalizes 
to the multi-operator case the framework of colligations (nodes) which has been 
extensively used in the study of a single nonselfadjoint (or non-unitary) operator, 
see, e.g., [3]. 

Let H be a Hilbert space (finite- or infinite-dimensional) and let E be a finite- 
dimensional Hilbert space. 





which implies the theorem. 



□ 
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Definition 3.2. An operator node is a collection 

C=(A,H,$,E,a) 

where A : H — > H is bounded linear operator, <E> is a bounded linear mapping from 
H to E with the adjoint mapping <f>* : E — > H , a is bounded selfadjoint operators 
in E, such that $*cr$ = I (A - A*). 

Definition 3.3. A commutative vessel is a collection 

V = (Ax , A 2 , H , $, E, (Ti , o- 2 , j m , 7 ou * ) 

where A\ , A 2 : H — » iJ are bounded linear commuting operators, ^ is a bounded 
linear mapping from H to E with the adjoint mapping <f>* : E — > <7i ; (72 > 7 m ; 
7°"* are bounded selfadjoint operators in E, such that 7™ = —j m ^ out = _^ OM * ; 
and 

z 

7 m $ = (7l$^-(7 2 $A^ 

7 °"*$ = a 1 ^A 2 -a 2 ^A 1 

7°"* = 7 m + i(<Ti$$V 2 ~ tr 2 ***<Ti) 

For the simplest case, when a system of operators consists of a single operator and 
the discriminant curve is a line, these problems were solved by Kravitsky in 1979, 
|13j . Let us consider an operator node C — (A, H,$, E, a) and three polynomials 
Po(x), pi{x) and p 2 (x) of degree n, such that po(A) is an invertible operator. We 
define A x = pi(A)p„ l (A) and A 2 = p 2 (A)p 1 (A) 

Theorem 3.4. A collection V = (A 1; A 2 , H, E', a x , cr 2 , 7™, j out ), where 

(7i = S(p ,Pl) ® cr 
(7 2 = B(p ,p 2 )®<7 

7™ = B( Pl ,p 2 )®(j 
7 out = 7 m + «((7i$'$'V 2 - (7 2 $'$'V 1 ) 

is a vessel and the the discriminant curve of this vessel is defined by the equation 
det m (B(p u p 2 ) + Xl B{p 2 , Po ) + x 2 B( Po , Pl )) = 

Explicit formulas for an image of a plane algebraic curve given by a determinantal 
representation under a rational transformation were obtained using generalization 
of Bezout matrices in ^H]. These results allowed to define properly the rational 
image of a system of two operators where the discriminant curve is a plane algebraic 
curve, [23j- 

3.2. Quadrature domains. The last two decades have witnessed a renewed in- 
terest and constant progress in the theory of quadrature domains. Questions such 
as the constructions and parametrization of quadrature domains with prescribed 
distribution u, the algebraic structure of the boundary or various functional ana- 
lytic characterizations of quadrature domains have been successfully investigated, 
see [221, |JE]- One of possible ways to find a polynomial that defines the boundary 
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of a quadrature domain is to consider such rational transformations of a complex 
plane that the image is a Riemann surface equipped with an involution. 

The domain Q C C is called a quadrature domain if there exists a distribution 
u with finite support in f2 such that 

fdA = u(f), 

n 

for every integrable analytic function / in il. To be more specific, there are points 
Xj £ £1 and constants jjk, < k < m(j) — 1, 1 < j < m such that 

m(j)—l 



Jo. rr; rr^ 



where m(j) is the multiplicity of the point Xj. The quadrature domain is an image 
of the unit disk D under polynomial transformation defined by a polynomial q(z) 
and the boundary of the quadrature domain is given by the equation A(z,"z) = 0. 
To find the polynomial A we consider a rational transformation of the complex 

plane: z -> {q{z),q{\)}. 

Theorem 3.5. If a quadrature domain is an image of the unit diskH> under poly- 
nomial transformation defined by a polynomial q(z) then there exist three polynomi- 
als po(z), pi(z) and p2(z) such that q{z) — ~t§ and q{ = ) — j^7§j- The boundary 
of Q is defined by the equation det(B(pi,p2) + zB{p2,Po) + zB(po,pi)) = 0. 

For the proof see . 

3.3. Braid monodromy. In 1937 Zariski laid down the foundations for the braid 
monodromy of curves in CP 2 . The braid monodromy is a homomorphism between 
the fundamental group of a punctured disk and the braid group. We give here a 
short description of the braid monodromy. 

Let us recall that the braid group B n is the A4CQ (mapping class group) of the 
n-puncture disk. We distinguish some important elements in the braid group B n 
which are called half-twists. Let D be a closed disk and let K = {kx, ■ ■ ■ ,k n } C 
D \ dD. Choose u £ dD. Let a, b be two points of K. We denote K a ^ = K\ {a, b}. 
Let (7 be a simple path in D\(dDUK a i>) connecting a with b. Choose a small regular 
neighborhood U of a and an orientation preserving diffeomorphism / : R 2 — > C such 
that /(cr) = [-1, 1], f(U) = {z£C | \z\ < 2}. 

Let a(x), < x be a real smooth monotone function such that: 

f 1, < x < | 
a(x) = < 



0, 2 < x 

Define a diffeomorphism h : C — » C as follows: for z = re %v £ C let h(z) = 

For the set {z £ C \ 2 < \z\}, h(z) = Id, and for the set {z £ C | \z\ < §}, h(z) 
is a rotation by 180° in the positive direction. 

Considering (/ o h o (we will compose from left to right) we get a diffeo- 

morphism of D which switches a and b and is the identity on D \ U. Thus it defines 
an element of B n [D, K\. 

The diffeomorphism (/ o h o / -1 )|.d defined above induces an automorphism on 
iri(D \K,u), that switches the position of two generators of tti(D \K,u). 
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Definition 3.6. Let H(a) be the braid defined by (/ o h o /~ 1 )\d- We call H{a) 
the positive half- twist defined by a. 

Let C be a real curve in C 2 of degree n. Denote by pr\ : C — » C and by 
pr2 : C — > C the projections to the first and second coordinate, defined in the 
obvious way. For x £ C we denote K(x) the projection of the points in C which lie 
with x as their first coordinate to the second coordinate (i.e., K(x) = pr2(pr 1 ~ 1 (a;))). 

Let N C C be the set N = N(C) = {x S C | \K(x)\ < n} = {xi, • • • ,x p }. We 
restrict ourselves only to the cases where N is finite. Take E to be a closed disc in 
C for which N C E \ dE. In addition take D to be a closed disc in C for which D 
contains all the points {K(x) | x 6 E}. That means that when restricted to E, we 
have C C E x D. 

With these definitions in hand we may define the braid monodromy of a projective 
curve: 

Definition 3.7. Let C be a projective curve of degree n in CP 2 , L be a generic 
line at infinity such that \L n C\ = n, and (x,y) is an affine coordinate system for 
C 2 = CP 2 \ L such that the projection of C to the first coordinate is generic. For 
E, D, N defined as above, let M 6 dE D R be the base point of ~K\(E \ N), and let 
a be an element of ~K\{E \ N). To a there are n lifts in C, each one of them begins 
and ends in the points of M x K{M). Projecting these lifts using pr2 : C — * C 
we get n paths in D which begin and end in the points of K(M). These induce a 
diffeomorphism of n\{D \ K{M)) which is the braid group B n as defined earlier. 
We call the homomorphism (p : -K\ (E \ N) — > B n the braid monodromy of C with 
respect to L, E x D,pri, and M . 

It is natural to ask questions about the connection between the rational trans- 
formation and the braid monodromy induced by its image. For example one can 
formulate the following questions: 

Question 3.8. Study the singular points of the image of a rational transforma- 
tion, and define conditions on po(x, y),p\{x, y),p2(x, y) which induce specific braid 
monodromy results. 

Question 3.9. Let C be a curve, let r be a rational transformation. Classify all 
braid monodromy results which may result. 

Question 3.10. Formulate necessary and sufficient conditions for a braid mon- 
odromy to be of rational curve. 

Question 3.11. Given a braid monodromy, which satisfies the sufficient condition 
above. Formulate a family of rational curves which will induce such braid mon- 
odromy. 

Question 3.12. Given two isomorphic rational curves. What can be said on their 
braid monodromies. 

Of course, these questions are more than wide, and at this point may not be 
completely answered. In |12| we established two results concerning degree 2 rational 
transformations, as follows: 

Let us consider the rational transformation 



(x,y) i ► (p (x,y),p 1 (x,y) 7 p 2 {x,y)). 
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In order to compute the local braid monodromy at the point (po{xo, yo), Pi(xo, yo), 
P2(x ,yo)), we assume that po(x ,y ) ^ 0. We define: n(x,y) = ^g^j , r 2 (x,y) = 

ll\x$) > and recursively 
D 1 (x)=r' 2 (x,0)- 7 r^ 

E 1 (y)=r> 2 (0,y)-^ 
E n (y)^D' n _ 1 (0,y)- T r^ ) 



Corollary 3.13. Let (po(xo,yo),pi{xo,yo),p2{xo,yQ)) be one of the intersec- 
tion points of the two conies at the image r(C). Let i be the minimal index for 
which Di(xo) ^ Ei(yo). Then, the multiplicity of the intersection point is i + 1, 
and thus the local braid monodromy at this intersection point is (i + 1) full twists 
of two strings. 

Theorem 3.14. |12| Let C be a curve which consists of two intersecting lines, and 
let r be a real rational transformation of degree 2 . Then, the braid monodromy of 
r(C) is completely defined by the number and multiplicity of it's real self intersection 
points. 

Theorem 13. 141 gives a full classification of the braid monodromy of the image of 
two intersecting lines under degree 2 rational transformations. 
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